Comparison of shear and dielectrics in highly viscous liquids 
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A relation between shear and dielectric spectra is derived for highly viscous liquids with a small 
rotational contribution Ae to the dielectric constant. It is valid if the shear fluctuations and the 
electric dipole fluctuations have the same spectrum. The comparison to literature data, taken under 
carefully controlled conditions to ensure samples from the same charge and the same temperature 
control in both measurements, indicates that the relation may be fulfilled or not depending on the 
substance. The connection to recent work on strong correlations is discussed. 
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A promising approach for the study of the flow process 
in highly viscous liquids is its comparisons^— in different 
techniques. One usually finds the dielectric absorption 
peak close to the heat capacity onei"— , but the shear 
modulus peak about half a decade higher in frequency. In 
a broad distribution of relaxation times, a modulus peak 
always appears at a higher frequency than a susceptibility 
(compliance) peak. The question is whether this is the 
reason here. 

The dielectric susceptibility e(a;) is difficult to invert 
because it has two contributions, the electronic polariz- 
ability and the molecular dipole orientation. It is better 
to invert the shear modulus G{lu). 

Let us start^ with the complex shear compliance 
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Here G is the infinite frequency shear modulus and 
T] is the viscosity. A third material constant hidden in 
this equation is the recoverable compliance Jg , the elas- 
tic compliance plus the integral over the retardation pro- 
cesses 



L{T)d Inr. 
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Eq. ID) makes a separation of two contributions to 
the compliance, the retardation spectrum and the vis- 
cosity. From our gradually growing understanding of the 
highly viscous liquid^'^, we know that both parts must 
come from thermally activated transitions between in- 
herent states, stable structures corresponding to minima 
of the potential energy. But the question is whether the 
retardation spectrum, which is only part of the response, 
can be used to calculate the normalized shear suscepti- 
bility from the fluctuation-dissipation theorem. 

In order to argue that this is indeed the case, consider 
an inherent state i of a small volume V in the liquid. V 
should be large enough to be characterized by the con- 
tinuum viscoelasticity, but small compared to the sample 
into which it is embedded. The state i has a well-deflned 
electric dipole moment Pi and a well-defined shear tensor. 
For simplicity, we consider only one of the five possible 



shear directions. Let be the shear strain of state i 
in this direction in the absence of external forces. The 
contribution Si of the state to the shear fluctuations of 
the volume V in the given direction depends on the time- 
dependent actual zero point eo{t) of the strain of the sur- 
roundings, which can take all possible values in a liquid. 
The actual strain of the volume V in the compromise of 
internal and external stresses is close to (e^ — eo{t))/2, 
depending on the Poisson ratio*. For our argument, one 
can take it to be exactly half of the difference. Then, 
the free energy Ei of the state i has a stress contribution 
GV{ei - eo)V4 and 

5. = /.(e.-eo(i))V4 (3) 
where the occupation probability fi is given by 

f^^^cM-E^/kBT), (4) 

Z being the partition function of the volume V for the 
given eo{t). 

In the case of a stationary flow cq ~ const, one has 



dE, _ 1 
dt 2 



GV{e, - eo{t))eo, 



(5) 



from which one can calculate the contribution eoi of state 
i to the flow 



GV 1 
2kBT Z 



Mei - eo(i))^eo. 
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The contribution of the state to the viscosity is propor- 
tional to its contribution to the shear fluctuations, with a 
proportionality constant which is the same for all states. 
This intimate relation between viscosity and fluctuations 
shows that the viscous flow is not a separate process. 

Thus the normalized shear susceptibility is deflned by 
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With J{uj) — 1/G{uj) and eq. ([T]), one can translate 
this to 
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Xs is normalized with respect to the comphance step; 
the function begins with the value 1 at low frequency 
and ends with at high frequency. To see this, note the 
limiting values of G(w) 



and 



and 



lim G' = G 



lim G" = J^r/^cj^ 
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(the last two are equs. (34) and (35) of chapter 3 of 
Ferry's book^). 

As long as the rotational contribution Ae = eiow — £high 
(where ciow is the low frequency limit and Chigh is the 
high frequency limit of the dielectric constant) is small 
compared to 1 , one can calculate the normalized complex 
dielectric susceptibility function Xe from the measured 
complex e(w) via 



Xe = 



(12) 



If the flow relaxation processes affect shear strain and 
electric dipole fluctuations in the same way, 



e" 



G 
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This equation will be used in the comparison to 
experiment^. 

The comparison is done for four samples with a 
weak Ae: TPE (triphenylethylene), squalane, 1,4- 
polybutadiene and DC704, a silicon oil used in diffusion 
pumpsi. 

One way to compare is to fit G(w), then calculate the 
expected normalized dielectric susceptibility from this fit 
and compare to the dielectric measurement. The fit is 
better than the data themselves, because small devia- 
tions in G(w) at low w tend to explode in the calculated 
retardation spectrum (see, for example, the increasing er- 
ror bars in Figs. 3 and 4 of Schroter and Donth's paper^). 
The fit also supplies the three parameters G, rj and 
needed for eq. (|13p . Ae is the only parameter which has 
to be adapted to the dielectric data. 

The procedure does not remove the accuracy problem 
at low frequency completely; the position of the resulting 
dielectric a-peak remains rather sensitive to the fit pa- 
rameters for G{uj). The physical reason for this is the vis- 
cosity, which overshadows the low frequency relaxations 
in a shear measurement (but not in a dielectric measure- 
ment). Here, G{uj) was modeled in terms of the asymme- 
try model^, introducing as an extra parameter a width 
w of the cutoff at the critical barrier Vc in order to get 
the best possible fit. 

In TPE, the first example, the relation is fulfilled 
within experimental error. This is seen in Fig. 1. Fig. 
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FIG. 1; Validity of the relation in triphenylethylene (TPE) 
(a) Fit of dynamical shear data* at 260 K (b) Comparison of 
the measured dielectric damping at the same temperature to 
the one calculated from the shear fit. 
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FIG. 2: Validity of the relation in squalane (a) Fit of dynam- 
ical shear data^ at 172 K (b) Comparison of the measured 
dielectric damping at 172 K and 190 K to the one calculated 
from the shear fit at the two temperatures, respectively. The 
arrows show the positions of the peak in G"{oj). 
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1 (a) shows the fit to G{uj), Fig. 1 (b) the comparison 
of the dielectric nreasurement to the prediction. Other 
temperatures yield the same result. 

The same good agreement is found in squalane, our sec- 
ond example. Fig. 2 (a) shows again the fit to G{u!) at 
172 K, Fig. 2 (b) the comparison between predicted and 
measured dielectric damping at two temperatures, 172 
and 190 K. This case is particularly impressive, because 
the shift between shear and dielectric peak increases from 
half a decade at the glass transition to one and a half at 
higher temperatures. Nevertheless, one finds good agree- 
ment at all temperatures, though the secondary peak is 
clearly weaker in dielectrics than in shear. 

The third case, DC704 in Fig. 3, does not show the per- 
fect agreement of the two previous cases. The calculated 
dielectric a-peak lies a bit to the right of the measured 
one. However, one also finds a slightly higher negative 
slope in the e" data, indicating a decrease of the ratio 
of the coupling constants of the elementary relaxations 
to dielectrics and shear with increasing frequency. If one 
corrects with a decrease of 6 % per frequency decade, one 
recovers a good fit. Again, the result does not depend on 
temperature. 

The last case, 1,4-polybutadiene (more precisely, a 
mixture of 80 % 1,4-polybutadiene and 20 % 1,2- 
polybutadiene) in Fig. 4 shows strong differences, with 
the calculated dielectric peak an order of magnitude 
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FIG. 3: Peak shift in DC704 (a) Fit of dynamical shear data^ 
at 219.5 K (b) Comparison of the measured dielectric damping 
at the same temperature to the one calculated from the shear 
fit (continuous line). The dashed line is obtained assuming 
a decrease of the dielectrics-shear coupling ratio by 6 % per 
frequency decade. 



slower than the observed one. Here, the deviations are 
much too pronounced to be repaired by the simple scal- 
ing which worked for DC704. However, this complete 
breakdown is in fact understandable. 1,4-polybutadiene 
consists of two kinds of monomers, cis and trans, with 
the cis-monomers responsible for the dielectric signal. A 
recent thorough simulation study of 1,4-polybutadiene 
finds that their response is decidedly faster than the one 
of the trans monomers. 

Our findings are not completely independent on our 
choice of a model^ for G{Ld); taking other models, one can 
even fit G(w) with a diverging recoverable compliance^. 
However, the results can be checked independently by 
an alternative procedure. One can fit the dielectric data 
first with a Cole-Davidson function (adding a Cole-Cole 
function for the secondary peak if necessary) . With this 
fit, one calculates a normalized dielectric susceptibility 
from eq. (fT2|). Then, one can check whether this suscep- 
tibility describes the G{uj) data correctly according to eq. 
([8]), taking not only G, r] and Jg as free parameters, but 
also the time constant of the Cole-Davidson function. If 
this agrees with the one in the dielectric fit within ex- 
perimental error (about a tenth of a decade) and if the 
fit is satisfactory, one again concludes that the shear and 
dielectric fluctuation spectra agree. In our four cases, 
this procedure corroborated the results from the G{uj)- 
fit, thus supporting the choice of a shear model with a 
finite recoverable compliance. 

Note that the identity of dielectric and shear retarda- 
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FIG. 4: Strong differences in polybutadiene (a) Fit of dy- 
namical shear data^ at 180 K (b) Comparison of the mea- 
sured dielectric damping at the same temperature to the one 
calculated from the shear ht. 
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tion spectra found in two of our four cases has a dif- 
ferent reason from the one of energy and density fluc- 
tuation spectra in strongly correlated liquidsiiii^ with 
a Prigogine-Defay-ratio close to one. In that case, the 
identity of the spectra is due to the strong correlation of 
the fluctuations. In our case, dielectric and shear fluc- 
tuations are not correlated at all, because a vector can- 
not be correlated with a tensor for symmetry reasons. 
Equal spectra do not necessarily imply strong correlation, 
though strong correlation implies equal spectra. In fact, 
the consideration on inherent states at the beginning of 
this paper shows clearly that the contributions of a given 
state to the shear strain and electric polarization fluc- 
tuations must be completely uncorrelated, because the 
first depends crucially on the external shear eo{t) and 
the second does not. 

The equality of the spectra must rather result from 
the decay mechanism. If the local strain disappears 
completely in a transition between two inherent states, 
the same must be true for the local dipole moment pi. 
In particular, the full decay of the shear strain should be 
accompanied by a full decay of the dielectric polarization 
(counterexamples are the normal modes in polymers or 
the monoalcoholsi^) . 

Summarizing the experimental and numerical evi- 



dence, one has two substances where eq. ([T3)) is valid at 
all measured temperatures within experimental accuracy 
and two other substances where one has good reasons 
to attribute the deviations to the frequency dependence 
of the coupling ratio. While a final judgment requires 
the investigation of more cases, the results strongly favor 
the conclusion that shear fiuctuation spectrum and di- 
electric fluctuation spectrum agree in simple substances 
and disagree if there is a complication like the one in 
1 ,4-polybutadiene. 

The finding provides for the first time a solid basis 
for the detailed comparison of the two most important 
techniques in the study of highly viscous liquids. If the 
spectra agree, the two techniques are complementary, the 
dielectric data providing the exact shape of the spectrum, 
in particular at the low frequency end which is not well 
seen in the shear measurement, and the shear data sup- 
plying the three material constants shear modulus, vis- 
cosity and recoverable compliance which together with 
the spectrum describe the flow. If the spectra do not 
agree, one can look for the microscopic reason. 

The paper profited a lot from intense discussions 
with the Roskilde group, in particular Niels Boie Olsen, 
Bo Jakobsen, Kristine Niss, Tage Christensen, Albena 
Nielsen and Tina Hecksher. 
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